Non-analytic behavior of 2D itinerant ferromagnets. 
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We consider an ordered ferromagnet in the vicinity of a T = transition into a paramagnet. We 
show that the free energy and the transverse and longitudinal static susceptibilities contain non- 
analyticities which destroy a continuous second-order transition. Depending on the parameters, the 
transition either becomes first-order, or occurs via an intermediate spiral phase. 

PACS numbers: 



Introduction. In recent years, there has been a strong 
interest to understand hidden features of a T = fer- 
romagnetic transition in itinerant fermionic systems. A 
much studied Hertz-Millis-Moriya (HMM) 0, 0, H </> 4 
model of a ferromagnetic quantum criticality predicts 
that the transition should be continuous in all dimen- 
sions D > 1, with mean-field exponents, like a classical 
transition in D > 4 = 1 + z, where z = 3 is the dynamical 
exponent. However, it was realized [J, la, [y, |7| that the 
difference between quantum and classical cases is more 
than the change of the effective dimension - in the quan- 
tum case, the 4 and higher-order terms in the expansion 
in the order parameter field contain singular dynamic 
parts, which reflect the fact that fermions give rise to 
long-range, dynamical interaction between collective spin 
excitations in itinerant fermionic systems. These extra 
dynamic terms have been analyzed on the paramagnetic 
side of the transition, both in D = 3 (Ref. J5i) and D = 2 
(Ref. [8j), and have been found to give rise to two com- 
peting effects: (i) the expansion of the free energy in the 
magnetic field is non-analytic, and the non-analytic term 
favors a preemptive first-order transition to a state with 
a finite magnetization, (ii) the static spin susceptibility 
x{q) is non- analytic in q and becomes negative at some 
qo, signaling another preemptive instability, this time to- 
wards a spiral phase. Which instability of a paramagnet 
comes first depends on the interplay between the prefac- 
tors for the analytic A 4 and q 2 terms, but in any event, 
a continuous second-order ferromagnetic transition point 
is internally unstable. 

In this paper, we consider what happens when the sys- 
tem approaches a ferromagnetic quantum-critical point 
(QCP) from the ordered state. We show that the free 
energy of a quantum ferromagnet is non-analytic in the 
order parameter field A, and the non- analytic term fa- 
vors a first-order transition into a paramagnet at some 
A c . We analyze the static spin susceptibility, which in 
the ferromagnetic phase has non-equal longitudinal and 
transverse components, X\\{l) an( i XJ-(<z)j and show that 
both are non-analytic in momentum q and both become 



negative at a finite q inside the ferromagnetic case, when 
A becomes smaller that some critical A c ± and A c ii . The 
negative x\\.±{l) implies the development of the spiral 
order, either along the direction of the magnetization, 
or in a transverse direction. We relate A c , A c ^ and 
A c ii with the parameters of the model and argue that 
if A c > A c _l,A c |i, the transition is first order, otherwise 
the ferromagnetic phase first becomes unstable against a 
spiral. 

Our analysis is based on the Eliashberg-type consid- 
eration near a ferromagnetic QCP. Such approximation 
has been justified in (y, |7], |8[ and we assume that it is 
valid. For dcfinitcncss, below we only consider the case 
D = 2, where the non-analyticities are stronger than in 
D = 3. 

The point of departure of our analysis is the spin- 
fermion model near a QCP. It describes fermions in- 
teracting with their own collective excitations in the 
spin channel, described by spin variables S = c\ o- a pcp. 
The model does not assume a long-range order a 'pri- 
ori, and is described by the Hamiltonian with three 
terms: H f = J2k,a=l,l e * c %a c l*;°" which describes 
low-energy fermions with the dispersion e^ rs Ef + 
vp{k — kp), H s — ^2„Xo ( < z)SqS_ q , which describes 
collective bosonic excitations with a bare static prop- 
agator Xo(<z)j and the spin-fermion interaction term 

Him = 9 ( 1 / Ar )Ek,q,a,a' C k,a Cr "^' S q C k+q,a'- Within 

this model, x(<7:^m) = Xo(<7)/(l - 2g 2 Xo(q)R(q, ^m)) 
where Xo(q) is the static susceptibility of free fermions, 
and II((7, £7 m ) is the polarization operator (here and be- 
low we set the Bohr magneton fis — 1)- Near a ferro- 
magnetic transition, x(l) = 2^(ej?)/((5 + (aq) 2 ), where 
v{ef) is the density of states per particle at the Fermi 
surface {v{ef) = m/(2n) for Eo(k) — k 2 /(2m)), 5 — 
1 — (2gv(eF)) 2 , and the length a is proportional to the 
radius of the interaction [8( . The dynamic fermionic self- 
energy and the dynamic part of H(q, Q m ) are computed 
self-consistently within the model. The fully renormal- 



ized 



spin susceptibility in the paramagnetic phase 



Xij(q,tt m ) =Sij 



2v(e F ) 
5 + {aq) 2 + n{q,n m )' 



(1) 



where II(g, J7 m ) = II(g, fl m ) — H(q, 0). For free fermions 
U(q,a m ) = \n m \/^/{v F q)2 + W m . 

The ferromagnetic transition occurs when g = 
\/(2v(e F )). At larger g, 5 becomes negative, and the 
system develops a ferromagnetic long-range order with 
the magnetization M = (S z ) = (7V T - N{). We as- 
sume that such long-range order does exist and search 
for preemptive instabilities upon approaching QCP from 
the ferromagnetic side. 

Mean-field analysis At the mean-held level, the in- 
teraction term reduces to g(S z ) ^7J fe a signa c^ a c^. a , and 
the fermionic propagator becomes 



G*\(k, u m ) = iu m - (e* - M) ± Aq, 



(2) 



where Ao = g(S z ) = g{N^ — N±), and jj, = ii(Aq) 
is the exact chemical potential (to order Aq, /j.(Aq) = 
iu(0) — i/'/(2i/)Aq, where v' is the derivative of the den- 
sity of states at the Fermi surface). Because Ao is finite, 
the longitudinal and transverse spin propagators become 
unequal already in the static limit. We have 



x zz (q) 

X xx (q) 



X z o z 



xV_ 

i-(<?x§ x ) 2 + M 2 ' 



(3) 



where 



Xo = - 
Xo x = -2 



duj m 

~2^~ 

dw, 



2tt 



v{e)de (G|(e,w m ) + G 2 (e,w m )) , 
v(e)dsG 1 (e,uj m )G l (e,u m ). (4) 



A simple calculation shows that x xx {fy = (JVf — 
N±)/A = 1/g, such that x xx (q) diverges at q — ► 0, in 
agreement with the Goldstone theorem. For longitudi- 
nal susceptibility, x zz ( a ) — '^ 1 '( £ f)/(S F + ( a q) 2 \ where 
5 F = -25 = KAl, and K = -(2/3)v"/v (Ref.9). The 
theory is only valid when K > 0, otherwise the transition 
is first order by trivial reasons. By order of magnitude, 
6f ~ (A /e F ) 2 . 

The dynamic terms U xx (q,fl m ) and H zz (q, f2 m ) also 
differ at Ao 7^ 0. Evaluating them using fermionic prop- 
agators from ([2]), we obtain 



n zz {q,£l m ) = 



n**( q ,in m ) = 



IO-. 



y/(v F q)* + ni 

I Or, 



V( v fQ) 2 + (fim + 2iA ) 



(5) 



The mean-field dynamic spin susceptibilities in the fer- 
romagnetic phase are then given by 



X zz (q,Vm) = 

X xx {q,itt m ) 



2v{e F ) 



5 F + (aq)2+Tl zz ( q ,n m )' 
2v{e F ) 



(aq) 2 +Il**{q,n r 



(6) 



Eliashberg theory Eqs. (J2J1B]) constitute the mean- 
field description of the ferromagnetic phase. Within this 
description, the transition is continuous, i.e., the ferro- 
magnetic phase is stable up to a point where Ao , 5 — 
0. The Eliashberg theory goes beyond this approxima- 
tion - it self-consistently takes into account w-dependent 
fermionic self-energy, but neglects k— dependent self- 
energy and vertex corrections. Vertex corrections gener- 
ally are not small if the interaction involve small momen- 
tum transfers, and are necessary to satisfy Ward iden- 
tities related to the conservation laws. However, the 
analysis of vertex corrections on the paramagnetic side 
have shown [7|, |8| that they can be rigorously neglected 
in the calculations of the non-analytic terms in the free 
energy and spin susceptibilities, if the interaction is suffi- 
ciently long-ranged such that ak F > 1, which we assume 
to hold [Hf • 

The calculations proceed in three steps, like in the 
paramagnetic phase. First, we obtain self-consistent one- 
loop expressions for the fermionic self-energy and dy- 
namic spin susceptibilities. Second, we use these one- 
loop expressions as inputs, obtain the free energy within 
Eliashberg theory, and show that it is non-analytic in 
A. Third, using the same inputs, we compute static spin 
susceptibilities at the two-loop level and find terms which 
are non-analytic in momentum. We argue that the non- 
analytic terms in the free energy favor a first-order tran- 
sition, while the non-analytic terms in the susceptibilities 
favor an intermediate spiral phase. 

Fermionic self-energy The one-loop fermionic self- 
energy in the ordered phase is given by 



S/(w„ 



A(£fc) LU m for LU m < LU /X 3 
LU ' W„{ for UJ m > LOq/X 3 



where ujq = 3^/3e F /(4:(ak F ) 4 ), e F — v F k F /2, and A(£fe) 
depends in non-singular way on the ratio of £& and Ao. 
For Ek = (J- 



A 



VF 



4ak F \aA ^J5 F + (2aA /v F ) 



(7) 



The self-energy is linear in uj m at the smallest frequen- 
cies, but crosses over to the quantum-critical, UJm be- 
havior at frequencies larger than uj /X 3 ~ Aq/(e f (ak F )). 
Such self-energy does not destroy the ferromagnetic or- 
der and preserves a Fermi surface, but it is larger than 
u m near QCP, and has non-Fermi liquid form in between 



ujq/X 3 <C wo and ojq. The non- Fermi liquid behavior in 
the ferromagnetic state has been discussed from a differ- 
ent point of view in Ref. ll| . 

The one-loop dynamic polarization operators 
nf x (<7, O m ) and n zz (g, fi TO ), re-evaluated with dressed 
fermions, are given by rather complex expressions. 
Like in previous studies [7], |8|, we found that, for the 
calculations of the non-analytic terms in the free energy 
and spin susceptibilities, we only need terms up to order 
1/q in the 1/q expansion of n(g, O m ). Such terms are 
not affected by vertex corrections [2|. We obtained 



nr(?,«) 

flffeo;) 



M 

M. 

v F q 



1 - 



(w + c a) S / (cj) + 2zA ) : ' 

2v 2 F q 2 
(w + c^E/^)) 5 



2«£g 



2 _2 



(8) 



1.20 



where c^ interpolates between c w= o = 1 and c u 
for ojq/\ 3 < lo < ujq. 

The free energy The free energy per particle for 
the ferromagnetic spin-fermion model in the Eliashbcrg 
approximation is given by 



S = S (A 



o;-t-=. 



2S* 



where S (A) = 2g& \-<l/4)5 F A 2 + (K/8)A 4 + ...] is 
analytic in A Minimizing S and expanding around 
the minimum, we obtain the equilibrium A = A = 
(Sf/K) 1 / 2 , and reproduce the mean-field expression for 
the static x zz (q ~ * 0)- Further, 



V 

2N 
V 
2N 



In 



dQ m d 2 q ,_ x**(o,o) 

(2tt) 3 
drj rri d 2 q' 



(27T) 



In 



X z *(<7,On) 

x xa (o,Q) 
x ra (g,fim)' 



(9) 



where x"(g, ft TO ) include one-loop polarization operators 
U- l i(q, Om)- Both S zz and S xx contain analytic contri- 
butions which renormalize constants in So (A). These 
renormalizations are small in l/(akp) (Ref. [8J]) and we 
neglect them. In addition, S 3131 contains the non-analytic 
term in A, which is our primary interest. Substituting 
\ xx with the polarization operator from (j8]) into Eqn. 
dH), integrating over momentum, and neglecting regular 
terms, we obtain at QCP 



y/2 A 7 /2 



2ttc 3 / 2 



1/2 



Z. 



(10) 



£ > w o 



where Z is the universal (cutoff independent) part of the 
integral 



Z = 2 



/ dyRe 
Jo 



tfl* -if logtf'*- if 



(11) 



The evaluation of the integral yields Z = — 87r\/2/35 « 
— 1.02. Combining non-analytic and analytic terms, we 
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FIG. 1: The free energy 5(A), Eq. JT2J, for different values 
of 6f, going from positive (6f)i to negative (<5_f)2,3,4- The 
first-order transition occurs at (Sf)3- 



obtain in the immediate vicinity of the QCP 



1 



-^A 2 + -A 4 
4 8 



A 7 /2 



£ 3 / 2 



(12) 



where E = dsp/iakp) 4 ^ 3 , and d rs 3.50. Apart for 
a small numerical difference in d, this expression coin- 
cides with the one obtained in [8j, where QCP was ap- 
proached from the paramagnetic side. There is, how- 
ever, an important distinction between paramagnetic and 
ferromagnetic phases away from QCP. In the paramag- 
netic phase, the A 7 / 2 dependence is replaced by A 3 when 
A/ef ~ S/(akp) 2 < S. In the ferromagnetic phase, 
S — ► Sf, which by itself scales as (A/ef) 2 . As a re- 
sult, the non-analytic A 7 / 2 dependence of S survives in 
the wide range away from a QCP, as long as A/sp < 7, 
where 7 = (ak F ) 2 /(Ke 2 f ). 

The negative A 7 / 2 term in the free energy shifts 
the equilibrium value A such that it remains finite 
even on the paramagnetic side of the transition, when 
Sf changes sign and becomes negative. For a generic 
S F , A is the solution of KAl - 7(A /E) 3 / 2 = 5 F - 
One can easily verify (see FigQ]) that the free energy 
(fi"2|) describes a first-order transition to a paramagnet 
at S F = -(A o /£0 3/2 < (i.e., 2<5 = -Sf > 0). 
The value of the equilibrium Ao at such transition is 
A c = 36/{K 2 E 3 ) « 0.847 2 e F . We also note that the 
stiffness - a prefactor for (A — Aq) 2 /2 in the free energy 



eff 



- changes from Sp to S 

and remains positive for all A 



eff 



S F + (7/4)(AoAE) 3 / 2 , 
> A c . At A = A c , 



(3/4)1*, 



Static spin susceptibilities We next show that 

the static spin susceptibilities x xx (q,0) and X zz {<lity 
also display non-analytic behavior, and that these non- 
analyticities compete with the one in the free energy and 
may give rise to pre-emptive spiral instabilities. The non- 
analytic term in %(g, 0) has been previously analyzed on 
the paramagnetic side [7|, |8| . We performed the calcula- 
tions in the ordered phase. 

The non-analytic behavior of x(g, 0) originates from 
non-analytic q— dependencies of II ZZ and H xx , which ac- 
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FIG. 2: The scaling function ^/yF xx (y). It has a maximum 
value of ~ 0.24 at y « 0.74 
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FIG. 3: The three critical A c as functions of the parameter 

j=(ak F ) 2 /(Ks F ). 



quire static parts at the two-loop order. 

The computational steps are similar to those in Ref. [7| 
and we refrain from discussing them. The non-analytic 
contributions IT; 2 and IG-p come from the processes in 
which fermions and spin fluctuations are vibrating near 
a fermionic mass shell and are far away from a bosonic 
mass shell [6j, |7| ( the same processes lead to fermionic 
self-energy E/(w m )). We found that the non-analyticity 
comes from the exchange processes involving transverse 
spin fluctuations. 

There are two types of non-analyticities in an ordered 
ferromagnet. First, there are corrections to static, uni- 
form II zz (0,0). They change S F into S F , which is the 
same as the stiffness obtained by expanding the free en- 
ergy. Second, there are non-analytic terms in the mo- 
mentum expansion of the static II zz (g,0) and H xx (q,0) 
For the latter, we found 



which is the solution of 



-a-q-'-k^, -F zz (2A /v F q), 



Il z 2 z (q,0) = -a 2 q 3/2 k 1 F /2 V zz < 

nf»(g, 0) = -a 2 q 3 / 2 k F /2 F xx (2A /v F q), (13) 



such that 



X zz (q,0) 



x xx (q,o) 



2v{e F ) 



2v(ef) 



1 



^- + g»-^ k ^ F "(^ 



q- - q", n, p 

1 



„2 _ „3/2I. 1 / 2 pii f 2A 

q q k f r y VFq 



(14) 



For A < v F q, F zz (0) = F xx {0) « 0.25, in agree- 
ment with [7|. For v F q <C Ao the two scaling func- 
tions differ. The scaling function F zz (y) remains close 
to F zz (0) as long as the argument y < e F /A . The 
function F xx (y) crosses over to F xx (y 3> 1) rs 0.15/^/y 
such that H xx (q, 0) becomes analytic, H xx (q, 0) ex q 2 . 
The function ^/yF xx (y) is plotted in Fig. 2 Analyzing 
Eq. (fl"4|) , we find that both susceptibilities become neg- 
ative at a finite q when Ao reduces below some critical 
value. The transverse susceptibility becomes negative at 
Aq = A Ci j_ = £F{^/yF(y)) 2 max = 0.06e F , whereas longi- 



e F 



1 - 



21 



7 



4 d 3 / 2 V A 



e F 



1/2 



= 0.02^7- 



Whether the transition is the first order or involves an 
intermediate spiral phase depends on which of A c , A c n 
and A C] j_ is the largest. All three critical A scale with 
e F , but they depend differently on the parameter 7. In 
FigJ3]we plotted the three critical A vs 7. We see that 
for 7 < 0.26, A c ,j^ > A Ci ||, A c , and the system first 
develops a transverse spiral order, while for 7 > 0.26, 
A c > A c m, A Cj _l, and the transition is first order. This 
qualitatively agrees with the analysis on the paramag- 
netic side [8|. 

In summary, we analyzed the non-analytic terms in the 
free energy and in static susceptibilities in an ordered 
itinerant ferromagnet. We found that, because of these 
non-analyticities, the transition to a paramagnet is either 
first order, or involves an intermediate spiral phase. 

We thank A. Andreev, P. Fulde, A. Green, A. Huxley, 
D. Maslov, A. Rosch, and Yu. Ovchinnikov for useful 
discussions. The research has been supported by NSF 
DMR 0604406 (A. Ch.) and in part by NSF PHY05- 
51164 (J.B.). J.B. would like to thank KITP, Santa Bar- 
bara for hospitality during the completion of this work. 



[1] 
[2] 

[3] 

[4] 

[5] 

[6] 

[7] 



tudinal susceptibility becomes negative at Aq 



^c,\\ 



J. Hertz, Phys. Rev. B 14, 1165 (1976). 
T. Moriya," Spin fluctuations in itinerant electron mag- 
netism", eds Springer- Verlag, 1985. 
A. J. Millis, Phys. Rev. B 48, 7183 (1993). 
D. Belitz, T. R. Kirkpatrick, and T. Vojta Rev. Mod. 
Phys. 77, 579 (2005) 

D. Belitz, T. R. Kirkpatrick, and T. Vojta, Phys. Re v. 
B 55, 9452 (1997). 

A. V. Chubukov, C. Pepin, and J. Rech, Phys. Rev. Lett. 
92, 147003 (2004). 

J. Rech, C. Pepin, and A. V. Chubukov, Phys. Rev. B 

74, 195126 (2006). 

D. L. Maslov, A. V. Chubukov, and R. Saha, Phys. Rev. 

B. 74, 220402(R) (2006). 



[9] A somewhat different result for K is obtained if one 
uses RPA for the Hubbard model. In this case, charge 
component of the interaction also contributes to K, and 
K = {v'/uf - (2/3)i/>. 
[10] An alternative way to justify Eliashberg theory even if 
akF < 1 is to extend the model to M > 1 spieces of 



fermions, which separately interact with a collective bo- 
son. Vertex corrections are then small in l/M. 
[11] V.P. Mineev, Phys. Rev. B72, 144418 (2005); D. Be- 
litz, T. R. Kirkpatrick, and A. Rosch, Phys. Rev. B 74, 
024409 (2006). 



